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A Submerged Singulor‘*.Uy Meuhod for CalcuXatlng 
Potential Flow Velocities a' Arbitrary Near-Field Points 

Brian Maskcw* 

Ams Heaearoh Center^ NASA^ Moffett Fietdj California 94035 

Summary 

A discrete singularity method has been developed for calculating the 
potential flow around two-dimensional airfoils. The objective was to calcu- 
late velocities at any arbitrary point in the flow field, Including points 
that approach the airfoil surface. That objective was achieved and is 
demonstrated here on a JouUowskl airfoil. The method used combined vortices 
and sources "submerged" a small distance below the airfoil surface and 
incorporated a near-field subvortex technique developed earlier. Wien a 
velocity calculation point approached the airfoil surface, the number of 
discrete singularities effectively increased (but only locally) to keep the 
point just outside the error region of the submerged singularity discretiza- 
tion. The method could be extended to three dimensions and should improve 
nonlinear methods, which calculate interference effects between multiple wings, 
and which include the effects of force-free trailing vortex sheets. The 
capability demonstrated here would extend the scope of such calculations to 
allow the close approach of wings and vortex sheets '(or vortices) . 

The author wishes to express his sincere appreciation to Mrs. Opal J. 

Lemmer for developing the plotting routines for this work; these routines 
reduced the analysis effort considerably. 

Index categories: Aircraft Aerodynamics; Subsonic Flow; Potential Flow. 

*NRC Research Associate; now Senior Research Scientist with Analytical 
Methods, Inc., 100 - ll6th S.E., Bellevue, Washington 98004. 
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Nomenclature 


A = length of vortex sheet represented by a discrete vortex (approximately 
equal to the distance between the two control points adjacent to the 
vortex) 

distance along the vortex sheet between the kth and the k. + 1 th 
basic vortices (Eq. (5)) 

0 = angle parameter for positioning the vortices, Eq, (3) 

A0 *= increment in 0 between two basic voitices 

60 = increment in 0 between two subvortices 

S = surface length (Eq. (3)) 

s = surface length measured to a basic vortex from the beginning of a 

V 

• region (i.e., from 0=0) 

R = position vector for basic vortices 

r = position vector for subvortices 

a = position vector of a calculation point relative to a vortex 

a = modulus of etc. 

n = unit normal vector at the airfoil surface 

^,k = cartesiaix unit vector system (Fig. 1) 

H = height of a velocity calculation point above the local vortex sheet 
_V = velocity vector 

= components of the vortex- induced velocity vector in the ^ and Jc 
W ) 

directions, respectively (Fig. 1) 
r = vortex strength 

D = piecewise-constant doublet strength, i.e., strength of opposing 
vortices at each end of a doublet panel (Fig. 3) 
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o « source strength 

C, = lift coefficient 

Cjj “ drag coefficient 

Cj^ » pitching moment coefficient about the origin 

NBS » number of basic singularities 

NSV “ number of subvortices used on a segment between two basic vortices 

FNS ^ factor on the number of subvortices (Eq. (1)) 

NRF “ near-field radius factor, applied to A (Subsection 3.3) 

LDF == submerged depth factor, applied to A 
SSF = subvortex strength factor (Eq. (4)) 

1.0 Introduction 

This is the second paper from a work program aimed at removing the close- 
approach problem associated with vortex- lattice methods. The first paper^ 
described a subvortex technique by which the near-field problem area of a 
discretized vortex sheet could be reduced to a small region, and showed an 
application in a free-vortex sheet rollup calculation. The present paper is 
concerned with potential flow pressure calculations at arbitrary surface 
points on thick, tv?o- dimensional airfoils. 

The objective of this study was to predict the velocity at any arbitrary 
point in the flowfield, especially at points that approached the discretized 
vortex sheet. ^ Not only would this enhance the versatility of vortex lattice 
methods,^”** but it would extend the force-free wake calculations®"^® toward 
close-approach situations involving multiple components and their wakes. 

Such a capability would particularly benefit the analysis of high-lift 
configurations and the calculation of other close interference effects between 
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wings and vortex sheets (or vortices) such ns occur in configurations with 
leading edge or tip-edge vortices. 

Although the present paper deals with the near-field problem in two- 
dimensional flow, the extension for three-dimensional methods (particularly 
for methods having a force-free wake) is a major consideration throughout. 

The development of the discrete singularity technique is described in Sec- 
tion 2, and a study of the effects of various parameters in the model is 
given in Section 3. The calculations are based on a cambered Joukowski 
airfoil, and are primarily concerned with the detailed pressure distribution, 
but, as a means of indicating overall accuracy, coefficients of drag — which 
should be zero here — lift and pitching moment also are evaluated by integrat- 
ing surface pressure effects. 

2.0 Development of the Method 
2.1 Basic Considerations 

Existing surface singularity methods, ^5-18 ,^i,ich are based mainly on 
piecewise constant singularity distributions on plane panels, have proved very 
powerful tools for predicting the potential flow characteristics of many 
configurations. However, in common with vortex lattice methods, their sur- 
face velocity calculations are essentially restricted to the control points. 
The present objective of calculating velocities at arbitrary points requires 
the near-field characteristics of a higher-order representation for both the 
singularity distribution and the airfoil geometry. Several higher-order 
surface singularity methods have been developed recently, but they do 
not necessarily have the capability for calculating velocities at arbitrary 
surface points. Furthermore, high-order representation can be cumbersome to 
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apply to three-dimensional high-litt configurations with multiple force-free 
wakes. Also, it can be wasteful in computing effort when calculating 
velocities at even a small distance from the singularity sheet unless care- 
fully controlled far-field models are Included. 

High-order representation adds little refinement to calculations in the 
far-field, yet many essentially far-field velocity calculations are performed 
when relaxing the trailing vortices in force-free wake applications.®’^ 

On the other hand, the simplest model for far-field calculations is based on 
vortex-lattice theory} indeed, some very versatile methods have been devel- 
oped, but velocity calculations close to the discretized vortex sheets 
are restricted to special lines of approach.^ 

2.2 Submerged Singularities 

One solution is to keep the singularity model simple and to place it 
inside the airfoil surface."^ The airfoil contour is then treated as a 
streamline of the flow. Several flow calculation methods include internal 
singularities, either on the chord line or on the cambec line. In the present 
work, internal singularities were developed such that they were placed on a 
"submerged" sheet closely related to the airfoil contour shape (Fig. 1). 
Obviously, the upper and lower parts of the sheet had a crossover upstream of 
the trailing edge, and so the affected parts were replaced by a single sheet 
extending from the crossover to the trailing edge along the mean line. 

The effort to calculate three-dimensional wake rollup makes it prefer- 
able to use discrete singularities, but thin airfoils might require many 

‘1‘This idea was suggested by Dr. V. J. Rossow in the Large-Scale Aerodynam- 
ics Branch at the NASA Ames Research Center. 



Bingularities to innintoln an accurate contour. For example, Fig. 2 shows 
pressure values calculated directly at arbitrary points on a cambered 
Joukowski airfoil which was represented by 41 submerged discrete vortices on 
a cosine spacing (described in Subsection 2.4). The vortex strengths were 
solved after specifying tangential flow at a set of control points on the 
airfoil surface.^ The control points and vortices are indicated in the lower 
part of Fig. 2. The Cj^ based on circulation was 1.4% in error. The 120 
pressure calculation points are not related to the vortex positions, and so 
the calculated pressure distribution shows large oscillations about the exact 
line (but it is much better than the distribution with the vortices on the 
surface). The submerged depth was 0.4A, whereas a depth of at least lA 
would be required to reduce the errors to an acceptable level. ^ But the 
required number of vortices might then be unacceptably large (bearing in mind 
the three-dimensional case) especially if the close relationship between the 
submerged sheet and the airfoil contour is to be maintained. The calculated 
results downstream of the crossover indicate that such a relationship, i.e., 
separate sheets as opposed to a single mean-line sheet, might be important 
for this model. (The submerged depth is examined in Subsection 3.2.) 

2.3 The Subvortex Technique 

The subvortex technique^ offers an attractive solution to this problem 
because it increases the number of vortices, but only where and when needed. 
Also, it is a logical model to use with the submerged singularity idea; for 
practical reasons, the subvortex technique gives a small region close to the 
singularity sheet where velocity errors are still appreciable.^ This region 
is now enclosed in the airfoil contour. Thus, calculation points approaching 
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the airfoil surface never "see" the holes in the discretlKation because, 
locally, the number of subvortices Increoses to keep the point just outside 
the error region. As the calculation point moves from the surface, however, 
the vortex model quickly rcv#?(’ts to the basic discretization. The expression 
in Ref. 1 that controls the number of subvorticos has been modified to improve 
the variation in NSV as the calculation point moves along the surface. The 
new expression is 

NSV = FNS { lnc.6.^v-pavt-0l [l + ,(n/^ °in/NRW‘^] | 

where FNS, an (even) input parameter in the computer program, allows the 
density of the subvortex system to be varied. The number of subvortlcea is 
kept even to ensure that the control points, as well as the basic vortex 
locations, are midpoints in the subvortex system.^ An upper limit is placed 
on the number of subvortices, based on a submerged depth factor (SDF) , i.e., 

(^integer-par t-of |^) . 


This prevents a runaway condition for the number of subvortices used near the 
trailing edge. 

The same induced velocity expression is used throughout, i.e.. 


V c= r 



( 2 ) 


where U = a /a^, and W = -a /a^. The vector a = (a i + a k) is the posi- 
tion vector of the calculation point relative to the vortex being considered 
(basic vortex or subvortex) and a^ = The near-field treatment, there- 
fore, is complicated only by the interpolation for the subvortex positions 
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and strengths; in the present work, linear intorpoiation was tried initially 
for both quantities to keep the model simple. 

l,t\ Vortex Positions 

In the earlier work, ^ the vortices were positioned with equal spacing, 
but in the present study, the liaitial vortices (before submerging) were 
positioned on the airfoil surfoce using equal angle increments in a cosine 
equation applied to distance along the contour, i.c., the dlstoncc along the 
contour to the kth vortex point is 

a «=» S (1 - cos 6^) > 2 ) 

\ 2 

where S is the length of fjntour associated with the interval 0 < 6 < r. 

In this spacing system, half-angles separate the initial vortex positions 
from the control points where the boundary condition of tangential flow is 
specified. This is an adaptation of Lan*s work;23 in keeps the singularity 
strength distribution more uniform when passing through "difficult" regions 
such as leading and trailing edges and flap hinge lines, With this point 
distribution, the first control point is located at the trailing edge, and so 
the Kutta condition is applied by specifying the flow direction there, e.g,, 
the direction along the mean line. 

From their Initial surface positions, the vortices are submerged along 
the local normal to the surface by a fraction of A, i.e. , SDFA. The sub- 
merged depth factor, SDF, is constant over the whole contour except near the 
trailing edge, where it automatically decreases along the single sheet 
(Fig. 1) . The control points remain on the airfoil contour except in the 
region very close to the trailing edge; here, corresponding upper and lower 
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control points arc combined and moved to the moan lino, Hence, the model 
adjacent to th<. trailing edge rosombloa a camber-lino model, and alleviates 
the Bcnaitivity of a surface vortlcity method to trailing edge shape. Because 
of this modeling, there arc more control points than unknown singularities, 
and so the equations are solved in a least-squares sense, using a NASA-Ames 
routine^'* based on the Householder method,^'® 

For the initial model, the subvortlces arc placed on straight segments 
joining the basic vortices (Fig, 3), They are positioned with equal incre- 
ments ill 0; for example, between the kth and (k + l)th basic vortices, 
the subvortex position vectors ares 

— — 1) H* i “ 1, 2, . • . , NSV (^i) 

wheres " (cos - cos 0Sj^)/(cos - cos 5 ~ 0.5)60, 

and 60 » A0/NSV. The subvortex strength factor (l.e. , subvortex strength/ 
basic vortex strength) is 

SSF ° e^g ~ n^)S sln(6B/2) 

^ ^k 

This distribution of subvortices approximates to a linear vorticlty variation 
on the straight segments joining the basic vortices. 

2.5 Equivalent Piecewise Constant Doublet Distribaclon 

In three-dimensions, quadrilateral vortices are convenient for modeling 
arbitrary geometry conflgurations.0»^»5 The present study, therefore, is 
based on the two-dimensional form of that model, viz., opposing vortex pairs 
(Fig. 3) which are equivalent to a piecewise uniform normal-doublet distribu- 
tion. Such a model, forming a closed surface, requires one doublet panel 
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ocrengdh to bo specified, otherwise the system is indeterminate. Accordingly, 
the upper panel adjacent to the crossover (Fig. 3) is specified to have zero 
strength. The boundary condition equation associated with the control point 
above the specified panel is still included in the system of equations, >"nlch 
is solved in the least-squares sense (Subsection 2.^0* The resultant vortex 
strengths are 


■■k ” \ " 1. 2 N 

where Dj^ are the doublet panel strengths, i.e., the strengths of the 
opposing vortex pairs. (Note that has been assumed zero.) 

2.6 Initial Results 

singularity model was used to calculate surface pressures on the 
cc;nt>«red Joukowski airfoil consldereJ in Subsection 2.2. Incidence was 10“, 
and the total number of basic vortices was 46 after submerging to a depth of 
O.IA. The initial vortex positions (before submerging) were identical to 
those in Subsection 2.2. The near-field radius factor (NRF) was 5.0, 
following the work of Maskew,^ and the subvortex parameter FNS in Eq. (1) 
was 2. The pressure distribution calculated at the same points as in Fig. 1 
is quite good (see Fig. 4) except at the leading edge and near the crossover 
between the upper and lower submerged sheets. Just upstream of the crossover, 
the corresponding upper and lower vortices have almost identical sets of 
influence coefficients, and this produces ill- conditioning of the equations. 

In the solution, therefore, corresponding upper and lower vortices form 
increasingly strong opposing pairs as the junction is approached, and so a 
strong local flow is induced approximately in the direction of the mean line. 
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Although this Increases the upper and lower surface calculated pressures 
(Fig. ^0» the resultant strength between the corresponding vortices gives a 
smooth chordwlse load distribution and an accurate lift and pitching inoraenw, 
e.g., the calculated results from pressure integration using the trapezoidal 
rule are = 1.706^f and Cg " -0.5369, compared with the exact values 
1.6973 and -0.5391, respectively. The integrated drag coefficient (which 
should be zero) is -0.0098, i.ti. , an error of -0.57% of Cj^. 

The computing time for this calculation is four times that for the basic 
case with no subvortex technique. A small part of this increase (16%) is 
attributable to the larger number of vortices in the present case, viz., 

46 cf. 41. The smaller submerged depth here results in fewer vortices being 
removed in the crossover region. The computing time could be reduced by 
storing the subvortex position vectors and strength factors; these quantities 
were recomputed each time in. the present program. 

2.7 Addition of Source Singularities 

One way of explaining the pressure deviations near the subsurface 
croBSOvei. is that corresponding upper and lower vortices are trying to 
provide a thickness effect at the surface (as well as the lifting effect) from 
a small base. Source singularities, which are more suitable for providing 
thickness effects, were therefore combined with the vortices to remove the 
problem by providing a more suitable basic onset flow for the vortices.’ The 
source strengths are evaluated from a local linearized solution based on the 
relative slope between the airfoil surface and camber line: 



The source onset flow is a unit flow along the direction of the local 

camber line to facilitate the treatment of cambered airfoils. (A uniform 
onset flow was found unsuitable in such cases.) The model increases the 
computation time by 8%, yet the sources receive the same "subvortex" treatment 
as the basic vortices, i.e., Eq. (2) becomes; 

V B fl’CUi + - U}i)l 

“ L 2ir J 

being the velocity contribution from a combined vortex/source. (In three 
dimensions the source treatment would be a more involved e'-tension of the 
vortex equation.) 

The sources, then, provide the vortices with a basic onset flow in which 
the thickness effects are approximately represented and this is particularly 
beneficial near the leading and trailing edge (Fig . 5) . The pressure 
distribution from this model shows marked improvement (compare Figs. 4 and 5), 
but there is still a tendency for the pressure distribution to oscillate, 
particularly near the leading edge. Integrating the pressure distribution 
yields the following force and moment coefficients: = 1.7040 (0.4% error), 

C,. = -0,0069 (0.4% of Ct), and C„ = -0.5377 (0.26% error). 

2.8 Increase in the Density of the Subvortex System 

To reduce the pressure o* jillation near the leading edge, the number of 
subvortices was doubled (FNS = 4 in Eq. (1)). Figure 6 shows the resulting 
pressure distribution. Examination of the values shows the whole distribution 
is improved to the point that even the minor undulations, which can just be 
detected in the upper surface distribution in Fig. 5, are removed; but, though 
improved, the leading-edge overshoot in pressure is still there. Doubling 


12 



Uhe factor FNS Increases computation time by 53%, while the integrated force 
and moment coefficients show only minor improvements in accuracy. Increasing 
the number of subvortices to FNS « 8 did not significantly improve the 
distribution in Fig. 6. 

2,9 Cur-ved Subvortex "Sheet” 

In a further attempt to remove the small leading-edge problem, a higher- 
order interpolation scheme (i.e., biquadratic) was applied for positioning 
the subvortices. Equation (4) was therefore replaced with 


1 “ Hj,* 1 ~ (1 “ 1 “ 


where 


n^(l - 

i2 

N— ■ M.| ; 

82 (rii>ni>n 2 ) 


g " (hi - ^ - Hj.) 


-AS 


Hi = 


k-1 


_ ASj^ + ASi^^+i 


AS,. 


TI2 = 


AS, 


h2 

as in Eq. (4) . 


(This is a cubic curve formed by a linear combination of two overlapping 
quadratic curves.) 

The modified technique was applied to the same Joukowski airfoil as 
before, with FNS = 2, The resulting pressure distribution, sho\m in Fig. 7, 
indicates a small improvement in the leading-edge region (compare with Fig. 5 
and examine the distances between pressure values near the leading edge) . 

The computation time increased by 17% over the linear- interpolation case. 

This increase is larger than needed in practice because the biquadratic 


13 



interpolation was applied throughout the airfoil, whereas it is clearly not 
required over a substantial part of the contour (compare Figs. 5 and 7). 

In a practical method, the higher-order interpolation would be applied only in 
regions of high curvature and onty when the maximum number of subvortices is 
being used in the segment. 

The higher-order geometry routine slightly Improved the accuracy for the 
integrated lift and pitching moment (the new errors being 0.34% and 0.24%, 
respectively), but the drag error was more than halved from -0.0069 to -0.0028 
(i.e., 0.4% of C, to 0.17% of ) as a result of the improved pressure 
distribution near the leading edge. Further improvement at the leading edge 
might be achieved using a higher-order strength distribution for the subvor- 
tices in the high curvature region. This would allow a closer representation 
of the extreme variation in vorticity that occurs at the leading edge in the 
present case (see also Subsection 3.1). Significantly, the pressure docs not 
oscillate when the large peak is removed from the pressure distribution, 
e.g., Fig. 8 shows the pressure distribution for the same airfoil at zero 

incidence (C error in this case was 0.04% of , i.e., 0.0002). 

D u 

2.10 Pressure Integration 

So far, the integrated values for C^, and have been based on 

the trapezoidal rule applied to 120 calculated pressure values. A number of 
cases were computed for the configuration in Subsection 2.9 (a ® 10°) varying 
the number of calculation points but keeping the basic vortex/source model the 
same. The resulting integrated and errors are shown in Fig. 9. 

Minor variations are indicated for numbers down to 50 and even 25 in the case 
of lift and pitching moment. With only 25 calculation points on the airfoil. 
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however, the integrated drag coefficient has started to diverge; apart from 
this, the drag error is held below 0.5% of and the lift and moment coef- 

ficient errors are below 0.5%. Small undulations occur because of the chang- 
ing relationship between the calculation paints and the small pressure 
oscillations seen in Fig. 7. Pressure integration for forces and monents was 
originally chosen as an overall guide to accuracy of the pressure distribution 
(particularly in the case of drag) . The alternative, based on vortex 
strengths, appears less sensitive to the present modifications; for the case 
in Fig. 7, the lift based on circulation is only 0.1% in error. 

2.11 Constraint Function 

The vortex pair/doublet strength distribution has a very smooth form 
(Fig. 10) because it represents integrated vorticity over the surface. The 
smoothness (when plotted against doublet subscript k) is helped by the use 
of the cosine spacing (Subsection 2.4), which concentrates the singularities 
in the region of large vorticity gradient. A further contribution to the 
smoothness of the distribution is provided by the source singularities 
(Fig. 10). The resulting distribution can be represented by fairly simple 
interpolation functions with a view to reducing the number of unknowns. This 
principle has been widely used in the past, particularly in linearized 
theories using global functions covering, for example, a full wing span or 
chord. In the present work, a biquadratic constraint function (similar to 
Eq. (6) but with n based on the doublet subscript k Instead of surface 
length) was briefly investigated to reduce the number of unknowns when the 
present method is extended to the three-dimensional form. The function was 
applied piecewise in a way similar to that in the constraint function 



technique being developed in the NASA Ames POTFAN program. The piecewise 
application (as opposed to a global one) should be more flexible in treating 
complicated configurations in three-dimensional flow. The full matrix of 
influence coefficients is still formed, but elements of the matrix are com- 
bined at selected pivotal positions in accordance with the constraint func- 
tions. After the pivotal values are solved, the complete set of doublet 
strengths is generated again using the constraint functions. 

Figure 11 shows the calculated pressure distribution for the Joukowski 
airfoil represented by 46 vortices with the number of unknowns halved. (The 
corresponding full solution is in Fig. 7.) The C. value decreased to 
1.6783 (-1% error), and the computation time increased 14%; this increase 
accounts for the matrix manipulation time since the saving in the solution 
time at the present low level of unknowns is insignificant in comparison with 
the overall run time. 

Further reductions in the number of unknowns caused substantial losses 
in the overall circulation and hence large errors in the pressure levels. 

The problem appears to be that the piecewise application of the constraint 
function decouples a large part of the doublet distribution from the Kutta 
condition control point, and although the general shape of the distribution is 
maintained (as in Fig. 10) the starting level at the trailing edge (which 
determines the circulation) falls as the number of unknowns is decreased. 
Development of this technique is continuing, and promises to remove the 
initial problem. The technique should lead to significant benefits for the 
three-dimensional method. 
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3.0 Parametric Study 




A parametric study was performed to examine the following factors in the 
method: (a) the number of basic singularities, NBS; (b) the submerged depth 

factor, SDF; and (c) the near-field radius factor, NRF. Throughout the 
parametric study, the status of the method was as follows: (a) sources were 

present and their strengths solved in accordance with Eq. (5)j (b) the 
biquadratic interpolation scheme was used for positioning the subvortices 
(Eq. (6)); (c) the constraint function routine was not used, i.e, , the full 
solution was obtained directly; and (d) the factor on the number of subvor- 
tices, FNS in Eq. (1), was 2. The calculations were performed for the same 
cambered Joukowski airfoil as in Section 2. Incidence was again 10“ and the 
pressures were calculated at the same 120 points. The base parameters, when 
they were not being varied, had the following values: NBS = 46, SDF = 0.1 

and NRF = 5.0. Exceptions occurred in two situations: first, when varying 

NBS, the submerged depth was held constant and so the factor SDF varied; 
second, when varying SDF, the number of basic singularities (NBS) varied 
slightly because of the changing length of the single sheet near the trailing 
edge. 


3.1 Effect of Number of Basic Singularities 

The number of basic singularities (i.e., combined vortex/sources) , NBS, 
was varied from 19 to 92 using a constant submerged depth. This depth 
corresponded to an SDF value of 0.1 in the base case with NBS = 46. Figure 12 
shows the pressure distribution for NBS = 19; using so few vortices would 
clearly be an advantage in three-dimensional applications. The calculated 
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pressure values are surprisingly good over moat of the airfoil, but an oscil- 
lation is present near the leading edge. This pressure oscillation decreases 
as NBS increases, and virtually disappears v?hen NflS “ 92. This result indi- 
cates that a closer representation of the conditions near the leading edge 
might improve the solution when using a small value for NBS; since the model 
already represents geometric curvature, and since we have already tried 
increasing the effective number of singularities (Subsection 2,8), then an 
Improvement might be obtained by using a higher-order distribution for the 
subvortex strengths near the leading edge (see also Subsection 2.9). 

Figure 13 shows the errors in the calculated values for the lift, drag 
and pitching moment coefficients for different values of NBS. These quanti- 
ties undulate slightly because of the numerical integration scheme, coupled 
with the changing relationships between the (fixed) calculation points and the 
(varying) vortex/source locations. Nevertheless, the error levels appear to 

be bounded: and C„ are well within 1 % error, and the C_ error is 

L M D 

below 0.5% of C , even with only 19 basic singularities. 

The computation time does not vary much for small values of NBS, 
e.g., the time with NBS = 19 is only 7% lower than the time with NBS = 46. 

The reason is that with a constant submerged depth, the required number of 
subvortices increases as the number of basic vortices decreases. The times 
for forming the matrix of influence coefficients and for solving the equations 
do vary with NBS, however, and although these form only a small part of the 
total time (i.e., less than 25%) for low NBS values, when NBS « 91 these 
contributions cause an increase of 64% in total time over that for the base 
case (NBS =* 46). 
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3.2 Effect of Submerged Depth 


The submerged depth factor, SDF, was varied from 0.025 to 0.^( on the 
basic case. Two of the resulting pressure distributions for SDF « 0.025 and 
0.4 are shown in Figs, 14(a) and 14(b), respectively. Figure 7 shows the 
corresponding distribution for SDF “ 0.1. With the small submerged depth 
(Fig. 14(a)), the leading-edge pressure oscillation is made worse — possibly 
because the "linear vorticity" distribution is closer to the surface — but 
the pressures in the traillng-edge region are even better than in the basic 
case with SDF = 0.1 (compare with Fig. 7). 

The larger submerged depth (Fig. 14(b)) gives a smoother pressure 
distribution near the leading-edge — but with a too-high value; at the trail- 
ing edge, however, the pressure distribution has collapsed, A value for SDF 
between 0.05 and 0.1 (for NBS “ 46) gives the best pressure distribution. 
Computing time should be taken into consideration, however. As the submerged 
depth decreases, the number of subvortices must increase, and so computing 
time increases rapidly, especially for SDF values below 0.1, e.g., the time 
for SDF •= 0.05 is 35% higher than that for SDF >=> 0.1. 

Figure 15 shows how the integrated force and moment coefficients vary 
with submerged depth. As SDF increases, and errors become rapidly 

worse — probably because of the deteriorating pressure distribution near the 
trailing edge. SDF does not appear to affect appreciably, 

3.3 Effect of Near-Field Radius 

The near-field radius factor, NRF, when multiplied by the A value of 
a basic vortex, defines a circle centered on chat vortex. ^Thenever a velocity 
calculation point comes inside the circle, then that basic vortex is modified 



by the subvortex technique. The near-flclcl radius was examined for a simple 
case,^ and on the basis of that, an NRF value of 5 was used for initial work 
here. To examine its affect in the airfoil application, NRF was varied from 
1 to 7. Figure 16 shows the pressure distribution with NRF 1 — which is 
clearly too small. The calculated pressure distribution at the arbitrary 
points improves as NRF increases, but there is little visual change in the 
distributions for NRF values above about 3. Computing time decreases rapidly 
as NRF is reduced; a value of 3 instead of 5 for NRF gives a time saving of 
30%. 

Figure 17 shows the effect of NRF on the force and moment errors from the 
pressure integration. They show excellent convergence characteristics as NRF 
increases, although appears to be converging towards an error of the 

order of 0.5%. The error in C, based on circulation, however, converges 

> L 

towards zero. 


<!|.0 Conclusions 

The combination of a near-field subvortex technique with a concept that 
places the singularities inside the airfoil has resulted in a method by which 
accurate pressures (and velocities) can be calculated directly (i.e., without 
interpolation) at any arbitrary point on the airfoil surface. The method is 
essentially a numerical integration procedure, but, by approaching it via the 
vortex lattice model, a useful set of rules and automatic procedures has been 
developed which makes the method accurate as well as efficient when moving 
from near- to far-field regions. The calculations were enhanced by combining 
sources with the vortices. 


20 



I 


1 


I 


i 


The rcsulca obtained so JJm,* Indicate that the number of h ^Ic slngulorl- 
tios used to represent a section should be of the order of ^iO u. 50. However> 
the results also suggest that the use of a hlgher~order strengtii variation for 
the Bubvortlces in regions of high pressure gradient might allow the number to 
be decreased — possibly as low as 20. 

Bearing in mind accuracy and computing effort, the optimum values for 
the submerged depth and for the near-field radius would appear to be of the 
order of O.IA and 3A, respectively. Computing time penalties quoted herein 
for the various parametric changes and developments are, in general, on the 
pessimistic side. There is considerable potential for improving the teat 
method to reduce computation effort, particularly in connection with the sub- 
vortex system. 

The method can be extended to three dimensions for application to vortex 
lattice based methods, and should then allow close-approach situations 
associated with multiple components and force-free wake calculations. 
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X CONTROL POINTS 

Pressures calculated at arbitrary points on a Joukowski airfoil at 
.deuce. Model: as in Pig, 5, but with number of subvortices doubled 
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Fig. 9 Effect of the number of calculation points on the errors in the inte- 
grated force and moment coefficients. Basic case as in Fig. 7. 
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rzg, ±± pressures calculated at arbitrar^r points on a Joukowski airfoil a 
10 incidence. Model: as for Fig. 7, but basic number of unknowns halved 
by using a constraint function. 
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Fig, 15 Effect of submerged depth factor, SDF, on the error 
grated force and moment coefficients. Basic case as in Fig. 
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■ 'g. 17 Effect of near-field radius factor, NRF, on the errors in the 
integrated force and moment coefficients. Basic case as in Fig. 7. 



